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Electronic states of isostructural single-component molecular conductors [M(tmdt)2] (M = Ni, Au, and Cu) are 
theoretically studied. By considering fragments of the molecular orbitals as the basis functions, we construct a multi- 
orbital model common for the three materials. The tight-binding parameters are estimated from first-principles band 
calculations, leading to a systematic view of their electronic structures. We find that the interplay between a p7r-type 
orbital (L) on each of the two ligands and a pda-type orbital (Ma) centered on the metal site plays a crucial role; their 
energy difference controls the electronic states near the Fermi energy. For the magnetic materials (M = Au and Cu), we 
take into account Coulomb interactions on different orbitals, i.e., consider the multi-orbital Hubbard model. Its ground- 
state properties are calculated within mean-field approximation where various types of magnetic structures with different 
orbital nature are found. An explanation for the experimental facts in [Cu(tmdt)2] is provided: The quasi-degeneracy of 
the two types of orbitals leads to a dual state where localized Ma spins appear while the L sites show a non-magnetic 
state due to dimerization. On the other hand, [Au(tmdt)2] locates in the subtle region of orbital mixing. We propose 
possible scenarios for its puzzling antiferromagnetic phase transition, involving the Ma orbital in contrast to previous 
discussions mostly concentrating on the L sector. 

KEYWORDS: molecular conductors, single-component molecular metal, multi-orbital system, Hubbard 
model, first-principles band calculation, magnetism, metal-insulator transition, Mott insulator 
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1. Introduction 

Molecular crystals composed of one kind of molecular 
species showing electrical conduction, i.e., single-component 
molecular conductors (SCMC), 1 ' 2 have been revealed to bear 
novel electronic states. In particular, their multi-orbital na- 
ture has been recognized, which is different from the situ- 
ation in conventional charge transfer salts (CTS); There, in 
most cases only one molecular orbital (MO) contributes to 
their electronic properties. 3,4 In fact, the involvement of dif- 
ferent MO is a consequence of the molecular design 5 in or- 
der to realize SCMC: to make the energy difference between 
the frontier MO small enough, so that their energy bands can 
overlap when inter-molecular transfer integrals become suf- 
ficiently large. Metal complex molecules of the form M(L)2 
(M = metal, L = ligand) are suitable for this purpose. Their 
frontier MO are approximately the bonding and anti-bonding 
combinations of the pir wave functions from the two ligands. 
Large ligands lead to an effectively small transfer integral be- 
tween them and results in a small energy difference. Such a 
situation indicates a two-MO system, which is realized also in 
some CTS as notably discussed in the M(dmit)2 based com- 
pounds. 6 ' 7 

The first SCMC where metallic conductivity was reported 1 
is [Ni(tmdt)2] (tmdt = trimethylene-tetrathiafulvalene- 
dithiolate), 8 whose resistivity decreases by cooling down to 
lowest temperatures (T). A direct evidence for its metallic 
feature was given by the observation of three-dimensional 



Fermi surfaces by de Haas-van Alphen oscillations, 9 whose 
results are consistent with first-principles band calcu- 
lations. 9,10 Near the Fermi energy, ep, there exist two 
overlapping bands from different £>7r-type MO and ep crosses 
the overlapping area. Electron and hole pockets appear, due 
to the existence of even number of electrons in the unit cell 
consisting of one Ni(tmdt)2 molecule. This is the success of 
the molecular design mentioned above. It shows Pauli para- 
magnetic susceptibility, 1 while isotropic magnetoresistance 
suggesting spin effect is observed 11 whose origin remains 
unclear. 

Since the discovery of [Ni(tmdt)2], many related com- 
pounds have been synthesized. Among them, an isostruc- 
tural analog but with odd number of electrons per unit cell, 
[Au(tmdt)2], 12 has been attracting interest. It shows an anti- 
ferromagnetic (AF) phase transition with the transition tem- 
perature Taf = HO K, 13,14 which is exceptionally high 
among molecular conductors. An intriguing point is that in 
the resistivity, showing metallic T-dependence down to low 
T as well, no anomaly around Taf is found. 15 Furthermore, 
analysis of an NMR measurement 14 suggests the magnetic 
moment in this AF state to be rather large, of the order of 
1 fiB • These features are distinct from the formation of a spin- 
density-wave state due to the nesting of Fermi surface, as fre- 
quently observed in CTS, where anomalies in transport prop- 
erties appear and typical values for the magnetic moment are 
one order of magnitude smaller, or even less. 4 Such a mag- 
netic solution attributed to the pir bands is actually stabilized 
in first-DrinciDles calculations 16 ' 17 as well as in a mean-field 
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(MF) study of an effective Hubbard model, 18 facing difficul- 
ties in explaining these experimental facts. 

Recently, another isostructural member [Cu(tmdt)2], hav- 
ing odd number of electrons per unit cell, was successfully 
synthesized. 19 It shows a semiconductive behavior in contrast 
to the two compounds above, and exhibits an AF phase tran- 
sition at T A f = 13 K, much lower than in [Au(tmdt) 2 ]. 19,20 
The T-dependence of magnetic susceptibility above Taf is 
ascribed to the behavior of the one-dimensional (ID) spin- 
1/2 Heisenberg model with AF exchange coupling of about 
150 meV, 19 which is consistent with 1 H-NMR nuclear spin- 
lattice relaxation rate (T-f 1 ) indicating ID spin dynamics. 20 
In this compound, in contrast with the discussions above, a 
pda-type MO centered at the metal site is suggested to lie 
close to the two ligandsp7r orbitals, and mix substantially. The 
charge transfer from the pir-MO results in a nearly half-filled 
pda-band. 21 Then, the magnetic properties of [Cu(tmdt)2] are 
attributed to localized spins appearing on the pda-MO. 19,20 

In fact, the possibility that more than the two pir-MO are in- 
volved in the electronic states of SCMC was first proposed for 
[Au(tmdt)2] based on first-principles band calculations: 16 ' 17 
The pdcr- and pir-MO mix slightly when forming the elec- 
tronic band structure, whereas the latter plays the major role 
near ep. However, more recently, it was inferred from exper- 
iments that the orbital energy difference between these MO 
is modified upon cooling by unusual structural variation, en- 
hancing the mixing. 22 

Such multi-MO characters in SCMC can be captured by 
the effective model approach based on the tight-binding ap- 
proximation, which has been successful and now widely used 
in CTS; 3 ' 4 ' 23 the observed de Haas-van Alphen oscillations in 
[Ni(tmdt)2] are consistent with the tight-binding picture. 9 In 
ref. 18, we proposed that the basis sets for the effective mod- 
els of [Ni(tmdt)2] and [Au(tmdt)2] can be taken as virtual or- 
bitals whose wave functions are parts of the relevant MO near 
ep, rather than the MO themselves. In this paper, we extend 
our theoretical approach to the newly synthesized [Cu(tmdt)2] 
and seek for a systematic view of the electronic states among 
the isostructural family of [M(tmdt) 2 ] (M = Ni, Au, and Cu). 

In § 2, we set up our effective multi-orbital Hubbard model 
and derive the tight-binding parameters by fitting to the re- 
sults of first-principles band calculations. By considering a 
common set of basis functions for the three materials, a sys- 
tematic view of the electronic states is achieved. Essentially, 
the transfer integrals providing the structures of each band are 
similar among the members, and the orbital mixing is mostly 
governed by the energy difference between the pit- and pda- 
type orbitals. 

Then, in § 3, by treating the on-site Coulomb interactions 
within MF approximation we investigate the ground state 
properties of models for [Au(tmdt)2] and [Cu(tmdt)2]. We 
will see that the orbital mixing brings about phase diagrams 
showing different magnetic states when the Coulomb interac- 
tions on the two types of orbitals are independently varied. 
In particular, a slight enhancement of mixing in [Au(tmdt)2] 
suggested by experiments 22 results in drastic changes from 
our previous results: 18 The involvement of the pdcr orbital is 
suggested. 

Sec. 4 is devoted to discussions, especially on the magnetic 
transitions observed in the two compounds. Our results are 



carries the ID S = 1/2 localized spins, interpreted as a multi- 
band Mott insulator. We discuss possible situations realized 
in [Au(tmdt)2], from the viewpoint of doped Mott insulating 
systems due to orbital mixing. A summary is given in § 5. 

2. Effective Model 

The wave functions that we chose as basis sets for the ef- 
fective model for SCMC in our previous work 18 are localized 
on some portions of the molecules. They can be considered as 
fragments of the MO, and then called as fragment MO (fMO) 
in refs. 24-26 which we follow in this paper as well. 

The original motivation to consider such decomposition of 
MO was the results of the first-principles calculations. 16 ' 17 
The spin-dependent calculation for [Au(tmdt)2] indicates a 
stable AF pattern where the spins align oppositely within 
each molecule. To understand this unusual situation, it was 
speculated that a ligands pir-pir transfer integral is larger be- 
tween adjacent molecules than that within a molecule, which 
is supported by the analysis based on the fMO approach. 18 
These features are consistent with the molecular design 5 men- 
tioned in § 1 and imply that the fMO approach gives an in- 
sightful picture. Such discussions have recently been elabo- 
rated within quantum chemistry calculations, 24-26 for (TTM- 
TTP)l3 and [Au(tmdt)2]. It is shown that the two approaches, 
the MO and fMO pictures, can be transformed from one to 
another. We note that the fMO approach here and the so- 
called fragment molecular orbital method applied to huge 
molecules 27 ' 28 share common concepts. 

In the following, all the first-principles calculations, includ- 
ing those for MO of isolated molecules, are performed using 
the computational code QMAS (Quantum MAterials Simula- 
tor) 29 based on the projector augmented- wave method 30 with 
the generalized gradient approximation. 31 See refs. 17 and 21 
for their details. 



Ni(tmdt) 2 Au(tmdt) 2 Cu(tmdt) 2 




Fig. 1. (Color online) Molecular orbitals for M(tmdt) 2 (M = Ni, Au, 
and Cu) molecules, from which the energy bands in their crystals near the 
Fermi energy are formed. 16 ' 17 ' 21 Non-magnetic wave functions and energy 
eigenvalues together with the electron occupation for isolated molecules are 
shown. The rough spatial extensions of the Mcr, L, Mtt orbitals explained in 
the text are indicated. 
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2. 1 Molecular orbitals and fragment model 

The band structures of isostructural SCMC near ep are 
composed of several MO with similar characters upon 
chemical modifications; 2 ' 16 ' 17 this applies to the family of 
[M(tmdt)2] including the new member M = Cu. 19 ' 21 The four 
MO which mostly contribute to the electronic bands near ep 
are shown in Fig. 1. 

They can be approximately reconstructed by three kinds 
of fMO which we call here Ma, L, and Mtt. The Ma and 
Mtt orbitals are the p-d mixed wave functions roughly being 
an anti-bonding combination of the metal site d xy and d xz 
orbitals, and the surrounding S 2p orbitals, respectively. 32 The 
relevant atomic d orbitals are 3d for M = Ni and Cu, and 5d 
for M = Au. The L orbital is the pn orbital which is similar 
to the HOMO of TTF molecule embedded in the ligands (see 
ref. 2); There are two of them in one molecule, LI and L2, 
one for each ligand; they are equivalent due to the inversion 
center at the metal site. 

In the isolated molecules, the Ma orbital does not mix with 
the others from their symmetry, then it is a MO itself, i.e., the 
pda -MO mentioned in § 1. As can be seen in Fig. 1, the other 
three MO can roughly be described as linear combinations of 
L and M-zr orbitals 18 as, </>li — </>L2 + ci</>mtt, 0li + </>L2> and 
0li — 0L2 — C2(/>mtt, respectively, where c\ and C2 are some 
coefficients 24 (we omit renormalization factors). 

In the fMO scheme, we consider these three kinds of or- 
bitals as a basis set composing the band structures, and then 
for the low-energy effective model. The two-MO case men- 
tioned in § 1 corresponds to the situation where only the L 
orbitals are considered. In ref. 18, we chose {L, Mtt} for 
[Ni(tmdt) 2 ] and {Ma, L} for [Au(tmdt) 2 ] to reproduce the 
first-principles band structures near ep. Here, all {Ma, L, 
Mtt} are taken into account as a common set, in order to pro- 
vide a systematic view among the compounds. 

Our model Hamiltonian including local Coulomb interac- 
tions reads: 

H=Ho J cH\ n t, (1) 

(l,m) s 

+ (Am. nf ° + A M7r nf*) , (2) 

i 

Hint = £ {U L {n^n\l + n^n\l) + U M o n™°nff 

i 

where Ho and Hint represent the one-particle part, determin- 
ing the band structure, and the on-site interaction, respec- 
tively. In the first term of eq. (2), ti m denotes the transfer in- 
tegrals between fMO, where the sum is taken for inter- fMO 
pairs (I, m) including intra- and inter-molecular ones, and q s 
(cj s ) denotes the annihilation (creation) operator for all kinds 
of orbitals with fMO site index / and spin s =t or |. In the 
second term of eq. (2), A Mc r and Amtt are the orbital energy 
of the Ma and Mtt orbitals, with respect to the L level. The 
sum here, as well as in eq. (3), is taken for the molecule in- 
dex i, where the number operators are n° is = c°J c° is and 
n° = n^+n^ with orbital index o = Ma, LI, L2, orM7r. The 
intra -orbital on-site Coulomb interactions for the three kinds 



of fMO are denoted as Ul, Uyia, and Umk, respectively. As 
for the inter-orbital on-site interaction, we only include 
between Ma and Mtt for simplicity, considering that these 
two orbitals share the spatial extent while they are separated 
from the L orbitals. 33 

2.2 Fitting to first-principles band calculations 

The tight-binding parameters are obtained by numerically 
fitting to the first-principles band structures for non-magnetic 
state. In Fig. 2, we show the bands near ep, together with the 
fitted tight-binding dispersions. The four bands originate from 
the four MO, or equivalently, the four fMO; The unit cell con- 
sists of one molecule. As for [Au(tmdt)2], calculations were 
performed for both room-T and low-T (9 K) structures de- 
termined experimentally, due to the indication of structural 
variation upon cooling, 22 as noted above. In the calculations 
for [Ni(tmdt)2] and [Cu(tmdt)2], the room-T structure param- 
eters are used. 

One can see that the top band in [Ni(tmdt)2] and the bot- 
tom band in [Au(tmdt)2] are separated from the others. This 
is the reason we previously used three-band fits (two kinds of 
fMO). 18 On the other hand, in [Cu(tmdt) 2 ] all four bands are 
overlapping, requiring a four-band fit for a reasonable agree- 
ment with the first-principles band structure. Note that the to- 
tal band widths of the four bands are about 1.7 eV (M = Ni) 
> 1.6 eV (Au)> 1.3 eV (Cu). 

The fitted tight-binding parameters together with the orbital 
occupancies per site calculated from Ho, i.e., (n°) = n (o = 
Ma, L, Mtt), are listed in Table I. The listed orbital energies, 
A^, and A^, are obtained by fitting the energy disper- 
sions of Ho - By noting that the first-principles band structures 
are obtained self-consistently including the Hartree contribu- 
tions within the interactions, we can make a correspondance 
between the fitted values and the orbital energies in eq. (1) 

QC 18,26,34 
dS, 

A^ a = A M a + + ^f^Tr ~ U h U h /2, (4) 

= Amtt + %l7rnM7r/2 + U^fiMa ~ U L fl L /2. (5) 

The transfer integrals ti m show more or less similar val- 
ues for all three members, 35 which is due to the fact that they 
are isostructural, as noticed in our previous work. 18 The Ma 
orbitals have large ti m only for the A bonds along the [100] 
direction: they show a ID structure. The L orbitals, on the 
other hand, possess a two-dimensional network, where the 
mter-molecular dimers are formed by the B bonds. Their net- 
work is schematically shown in Figs. 3(b) and (c); in the unit 
cell, the two L sites from different molecules form the L1-L2 
dimers. The degree of dimerization, represented by the intra- 
dimer transfer integral, namely that along the L1-L2 B bond, 
is largest in [Cu(tmdt)2]. This dimerized structure resembles 
the situation commonly seen in typical CTS, e.g., in TM2X 
and polytypes of ET2X such as k and /3-type. 3,4 The Mtt or- 
bitals, in contrast, do not have large ti m between them, but 
bridge the L layers mainly along the intra-molecular bonds; 
They are appreciable (0.1 - 0.2 eV) and therefore these two 
orbitals mix with each other. The main differences between 
the three compounds are in the orbital energies, which we will 
discuss in the next subsection. 
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(a) [Ni(tmdt) 2 ] 
1.0 




-0.5 



-1.0 



r YT zs 
(c) [Au(tmdt) 2 ] (9 K) 



UT 




(b) [Au(tmdt) 2 ] (R.T.) 
0.5 




Yr zs 



Fig. 2. (Color online) The tight-binding dispersions (lines) obtained by fitting to first-principles band structure (dots) near the Fermi energy ep (dotted lines) 
of (a) [Ni(tmdt)2], [Au(tmdt)2] for (b) room temperature and (c) 9 K structures, and (d) [Cu(tmdt)2]. The dispersions are drawn between T-point, and X(7r,0,0), 
Y(0,7r,0), Z(0,0,7r), S(7t,7t,0), U(7t,0,7t), T(0,7t,7t), and R(7r,7r,7r)-points, and first-principles calculations data used for the fitting along them are plotted. 



Table I. Tight-binding parameters, ti m , and orbital energies from the L level, Ajjj^, and Ajjj^, obtained by fitting to first-principles band calculations. For 
example, bond: A, site pair: Ll-Mcr, denotes £/ m along the A ([1,0,0]) bond between LI and Mcr sites belonging to the molecule at the origin (0, 0, 0) and that 
at (a, 0, 0), respectively; As for notations for the inter-molecular bonds, see Fig. 3(a). ti m with large absolute values are shown by bold characters. Orbital 
occupancies per site, n (o = Mcr, L, Mtt) are also listed. 



bond 


site pair 


[Ni(tmdt) 2 ] 


[Au(tmdt) 2 ] 


[Au(tmdt) 2 ] (9 K) 


[Cu(tmdt) 2 ] 


intra-mol. 
A 
B 

Q 
P 


L1-L2 
L1-L1,L2-L2 
L1-L2 
L1-L2 
L1-L2 


-61 meV 
-86 
221 
131 

33 


-24 meV 
-99 
207 
109 

41 


-27 meV 
-98 
222 
129 

46 


-35 meV 
-90 
250 
134 

37 


A 
B 
C 
R 


Ll-Mcr, Mcr-L2 
Ll-Mcr, Mcr-L2 
Ll-Mcr, Mcr-L2 
Ll-Mcr, M(7-L2 


37 
27 
-50 
21 


33 
27 
-25 
22 


41 

33 
-29 
27 


14 
29 
-24 
14 


intra-mol. 
B 
C 

Q 


LI-Mtt, Mtt-L2 
LI-Mtt, Mtt-L2 
LI-Mtt, Mtt-L2 
LI-Mtt, Mtt-L2 


-200 



-22 
5 


-138 

30 
-5 
15 


-127 

35 
-22 
14 


-188 

20 
-7 
11 


A 
A 
B 
A 


Mcr-Mcr 
Mtt-Mtt 
Mtt-Mtt 
Mtt-Mct 


82 

-13 

24 
53 


95 

-23 
37 
-34 


98 

-30 
44 
-14 


96 

-19 
50 
-1 


orbital energies 


A u 

Mcr 

A 




879 meV 
-193 


522 meV 
-612 


476 meV 
-636 


181 meV 
-377 


orbital occupancies 


flh 
UMtt 




0.02 
1.16 
1.66 


0.03 
1.50 
1.97 


0.09 
1.47 
1.97 


0.86 
1.15 
1.85 



2.3 Systematic view of electronic structures 

The orbital energies A^ and A^ together with the trans- 
fer integrals determining the band structures, and the electron 
numbers, 4 (M = Ni) or 5 (M = Au, Cu) in the four orbitals, 



lead us to a systematic view for this family, as shown in Fig. 4. 

The Mir level is low in energy but mixes with the L bands 
especially in [Ni(tmdt)2] and [Cu(tmdt)2]; umtt deviates from 
2, as seen in Table I. However, as far as the main characteris- 
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[Ni(tmdt) 2 ] [Au(tmdt) 2 ] [Cu(tmdt) 2 ] 




Fig. 4. (Color) Energy diagram of orbital energies based on Table I, with respect to the Fermi energy ep, together with the approximate electron occupancies: 
4 electrons for M = Ni and 5 electrons for M = Au and Cu in total, in the unit cell of one molecule having 4 orbitals. The degenerated L level (broken red line) 
splits due to the inter-molecular dimerization; Aedimer represents the dimerization gap. Schematic density of states (DOS) are also shown for each orbital: Mo- 
has one-dimensional electronic structure, and the dimerization in L is respresented by the two hump structure. Note that the orbital mixing is not considered 
here. (See refs. 16, 17,21 for the actual DOS by first-principles.) 




Fig. 3. (Color online) [M(tmdt)2] viewed along the molecular long axis 
(a) and schematic views of the fMO model with lattice sites represented as 
circles [(b) and (c)] . The unit cell containing one molecule with four orbitals 
is shown as gray lines, and notations for the inter-molecular bonds are in- 
dicated whose correspondences are: A [100], B [111], C [101], P [211], Q 
[001], R [01 1]. The two-dimensional network of L sites is shown in (b), per- 
pendicular to the plane, and in (c); the sites connected by the largest transfer 
integral along the B bonds form dimers. See also Fig. 1 of ref. 18. 



tics of the electronic states near ep is concerned, it does not 
play an important role. Then, the value of is the essen- 
tial difference when comparing the three members, control- 
ling the orbital mixing between the Ma and L orbitals. It be- 
comes monotonically small as M = Ni Au (room-T) — >• 
Au (low-T) — >• Cu. This gives rise to a crucial difference in 
the magnetic states of Au and Cu systems as will be seen in 
Sec. 3. 



The Ma level in [Ni(tmdt)2] is high and approximately 2 
electrons enters the L level. The existence of two L sites in the 
unit cell which show dimerization as discussed in the previous 
subsection results in a band splitting, but not large enough 
to make a direct band gap. Then electron and hole pockets 
appear and compose the Fermi surface; This is how the first 
SCMC with metallic character was realized. 

Among the parameters for [Au(tmdt)2] comparing the 
room-T and low-T structures, the low-T data show slightly 
larger transfer integrals than room-T especially in the L-L 
pairs; this is naturally expected from the thermal contraction. 

is reduced by lowering T by about 0.05 eV, which is 
consistent with the previous MO calculation. 22 As a result, 
a fraction of the electrons in the L sites, nearly 3/4-filled in 
the room-T parameters, are transfered to the Ma orbital in 
the low-T parameters: 1%^, n^} = {0.03, 1.5} (room-T) —> 
{0.09, 1.47} (low-T). 

As for [Cu(tmdt)2], Ma mixes more with L due to the re- 
duction of further, about 0.30 eV smaller than the low- 
T parameter for [Au(tmdt)2]. Especially, near ep the con- 
tribution from Ma is appreciable, 21 even though the orbital 
level scheme shows positive A^. This is because the L or- 
bitals have larger ti m with a two-dimensional character while 
the Ma band is ID, therefore the former show wider bands. 
The orbital occupancies are {fiMa, ^l} = {0.86, 1.15}, which 
are substantially varied from the cases above. These features 
are close to the situation in the MO scheme in Fig. 1: In 
Cu(tmdt)2, the Ma (pda) orbital is occupied with nearly one 
electron. 

Let us comment on the relation between our level 
scheme and the nominal charges given as Ni 2+ [(tmdt)~]2, 
Au 3+ [(tmdt) 1 - 5 "] 2 , and Cu 2+ [(tmdt)~] 2 ; the correspond- 
ing occupations of metal d levels are Ni:(3d) 8 , Au:(5d) 8 , 
and Cu:(3cT) 9 , respectively. When one hypothesize the Ma 
and Mtt levels as the 'cT levels of the metal atoms, and 
the L orbitals for the full charge of the ligands (namely, 
omit the p-d mixing which is actually large), these nomi- 
nal charges correspond to, (M7r) 2 [(L) 1 ] 2 , (M7r) 2 [(L) 15 ]2, and 
(M7r) 2 (Ma) 1 [(L) 1 ]2, respectively; namely, a full-filled Mir 
for all compounds, and, (i) in M = Ni, a 1/2-filled L band, 
(ii) M = Au, a 3/4-filled L band, and (iii) M = Cu, a 1/2- 
filled Ma band and a 1/2-filled L band. These are close to the 
situations in Fig. 4. 
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(b) L-AF q mag [L]=(3t,0,0) 




Ma/Mjt 



(c) L&Ma-AF q mag [L]=q mag [Ma]=(^,0) 




( d ) [Au(tmdt) 2 ] 

a.b. 4 Ma 



,4f 



[Cu(tmdt) 2 ] 



Ma 



b. 



Fig. 5. (Color) Schematic views for (a) unit cell containing 4 orbitals with 
the red bar indicating the inter-molecular B bond in Fig. 3, and (b) and (c) rep- 
resentative antiferromagnetic (AF) solutions projected on the ab plane [see 
Fig. 3(c)]: (b) L-AF: Spin moments appear only on L sites [the dimer-AF 
pattern with spin ordering vector g m ag[L] = (tt, 0, 0)] seen for parameters 
of [Au(tmdt)2]; 16 ' 18 (c) L&Mcr-AF: Spins align in a staggered AF manner 
among the L and Ma sites in the case of [Cu(tmdt)2] with the common or- 
dering vector g m ag[L] = g m ag[Mcr] = (tt, 7r, 0). In (d), schematic represen- 
tations of electron occupancies realized in parameters for [Au(tmdt)2] and 
[Cu(tmdt)2] are shown, respectively. and 'a. b! represent the bonding 
and antibonding orbitals for L dimer. 



3. Mean-Field Calculation 

As mentioned in § 1, [Au(tmdt)2] and [Cu(tmdt)2] show 
phase transitions to magnetically ordered states. Here we 
study their ground-state (T = 0) magnetic states where the in- 
teraction terms in Hint [eq. (3)] are treated within MF approx- 
imation as n° i8 n? 8 , Such 
a MF treatment is suitable in seeking for possible different 
states, as in our model here with multiple degrees of freedom. 
In the calculations, we consider a supercell of 2 a x 2b x 2c, 
which includes {Ma, L (LI and L2), Mir} x 8 orbitals [see 
Fig 5(a)]. 

We set the values of tim to the fitted results in Table I, while 
Amct and Amtt are adjusted by the condition in eqs. (4) and 
(5) when the interaction parameters are varied. Then the band 
dispersions are unchanged within the paramagnetic metallic 
(PM) solution. Namely, the orbital occupations in the PM 
state, (n°^) = (n°^) = (n°)/2, are fixed at the values listed 
in Table I as, (n°) = n . As for the magnetic solutions, we 
relax the condition of fixed occupation and searched for self- 
ronsistent solutions of lowest enercrv in an unrestricted man- 



ner within the periodicity of the supercell. 

In the following, we set Uyia = Umtt = U M = Um for 
simplicity 36 and vary Um and C/l independently. In the whole 
parameter range we sought, the Mir orbital has negligible spin 
moment, i.e., it is magnetically inactive. Then the parameters 
Um and C/l control the correlation effect on the Ma and L or- 
bitals, respectively. A speculation we can make is the relation 
Um > Ul, considering the d contribution to the Ma orbital as 
well as its smaller spacial extent than the L orbital. 

Before presenting the results, we remark about two typical 
AF solutions. The first is stabilized in the case of [Au(tmdt)2] 
when C/l is increased, whose nature was discussed in refs. 16 
and 18. As shown in Fig. 5(b), spin moments appear only on 
the L sites therefore it is represented as L-AF. The spins are 
parallel within the L1-L2 dimers connected by the B bonds 
[see Fig. 3(c)], and anti-parallel between dimers along the 
inter-dimer bonds denoted as A, P, and Q in Fig. 3 (a); Their or- 
dering vector is g m ag[L] = (7r, 0, 0). This pattern corresponds 
to the dimer-AF spin order frequently appearing in 1/4-filled 
CTS under dimerization. 3 When the on-site Coulomb repul- 
sion is sufficiently large, their staggered pattern can open a 
gap at ep when the system is 1/4-filled in terms of either elec- 
trons or holes with each dimer carrying effective S = 1/2. 
Such a case is considered as the dimer-Mott insulating state. 
In the case of the room-T parameters for [Au(tmdt)2], the 
L bands are 3/4-filled, without pronounced mixing with the 
other orbitals. The dimers are formed by inter-molecular fMO 
sites and the spins on the two L sites in a molecule become 
anti-parallel. 16 The results of the first-principles calculation 16 
corresponds to relatively small C/l case and the system re- 
mains metallic even in the L-AF state. 

On the other hand, in the case of [Cu(tmdt)2], a stable AF 
pattern shows spin moments on both L and Ma sites when 
both Um and C/l are large enough, then denoted as L&Ma- 
AF. As shown in Fig. 5(c), the spins on the L sites show a 
staggered AF state where those connected with A, B, and Q 
(see Fig. 3) bonds are anti-parallel; ti m are large along these 
L-L bonds as shown in Table I. The spins on the Ma sites 
are also staggered, with the common spin ordering vector 
#mag[L] = gmagpVla] = (jr, 7r, 0). This can open an insulating 
gap when both L and Ma bands becomes nearly 1/2-filled; 
both orbitals provide effective 5 = 1/2 and the AF pattern 
corresponds to the Neel state configuration. Now this can be 
considered as a multi-band Mott insulator, due to the quasi- 
degeneracy of the two orbitals. 

The two contrasting situations for large interactions are 
summarized in Fig. 5(d), where schematic representations of 
electron occupancies in the two cases are shown. In the case 
of [Au(tmdt)2], the localized spins appear on the anti-bonding 
orbital of the dimerized L sites. In contrast, in the case of 
[Cu(tmdt)2] the spins appear on each of the L sites as well as 
the Ma sites as a result of orbital mixing. 

3.1 [Au(tmdt) 2 ] 

In Fig. 6, we show ground- state phase diagrams for 
three different parameter sets corresponding to [Au(tmdt)2]. 
Figs. 6(a) and (b) are those with the fitted results for the 
room-T and low-T structures, respectively, in Table I. Be- 
sides them, to see the effect of reduction of A Ma more ex- 
plicitly, we artificially decrease it from its room-T value as 
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rameters unchanged: This is shown in Fig. 6(c). The distinct 
states are indicated in the phase diagrams by different sym- 
bols summarized in Fig. 6(d), together with their magnetic 
ordering vectors. As seen there, the AF order on the L sites 
always has the pattern in Fig. 5(b) with gWg[L] = (7r, 0, 0). 
By the reduction of A^, we indeed find a crucial difference: 
The dimer-type AF insulating (AFI) state within the L orbital 
discussed above [filled symbol in Fig. 6(a)] stabilized in the 
wide range of parameters in the room-T case is not seen, and 
different AF metallic (AFM) states appear in Figs. 6(b) and 
(c), due to the mixing between the L and Ma orbitals. 

In Fig. 7, we show the /7 L -dependence of site occupation 
number for the three kinds of orbitals together with the spin 
moments on the L and Ma orbitals. Figs. 7 (a)-(c) correspond 
to the cases in Figs. 6(a)-(c), respectively. The right and left 
panels in Figs. 7 (b) and (c) are data along different traces 
in the (/7m, E/l) plane, i.e., for Um = Ul and Um = 2Z7l, 
respectively. 

3.1.1 Room-temperature structure 

The phase diagram in Fig. 6(a) shows no dependence on 
Um', this is because the Ma orbital is always nearly unoccu- 
pied, as seen in Fig. 7(a). As is increased, the system varies 
as PM — ^ AFM — >> AFI states. Since only the L sites possess 
magnetic moments the magnetic phases are represented as L- 
AFM and L-AFI in Fig. 7(a). 

These results are almost identical to the results in ref. 18, 
where the MF calculations were carried out for the three- 
band model based on {Ma, L} and assuming Umo- = U^. 
This is consistent with the fact that present results show an 
almost fully occupied Mir orbital playing no role, and no de- 
pendence on Um- As discussed there, the properties are gov- 
erned by the L orbitals. Their occupation number of 1.5 to- 
gether with the rather strong dimerization, are consistent with 
the fact that our results are analogous to the MF calculations 
on two-dimensional 1/4-filled Hubbard models with dimer- 
ization, e.g., as firstly performed on the model for ft-ET2X. 37 
The existence of the AFM phase in between the AFI and PM 
phases is due to the imperfect nesting property of the Fermi 
surface, where small magnetic moments cannot produce a 
band splitting large enough to open up a gap on the whole 
Fermi surface. As far as the system is in the AFM phase, the 
spin moment on each site is less than 0.2/i B ; this is also the 
same as in the three-band model.. 18 

3.1.2 Low -temperature structure and reduced A^ CT 

The phase diagrams, Fig. 6(b) for the fitted results for the 
low-T structure, and Fig. 6(c) for the room-T values but 
with reduced A Mcr , share common features in their overall 
structure, and are very different from the room-T case. The 
spin moments appear on the Ma orbitals when Um is en- 
larged, while insulating states are not stabilized, at least for 
(Ul,Um) < 0.8 eV; the whole phase diagrams show metallic 
states. The phase diagrams are devided into four regions: PM, 
L-AFM, magnetic metallic states with moments on the Ma 
sites [either antiferromagnetic (Ma- AFM) or ferromagnetic 
(Ma-FM)], and those with moments on both orbitals [L&Ma- 
AFM or L&Ma-FM]. These are common for Figs. 6(b) and 
(c), therefore the main variation from room-T to low-T struc- 
ture can be captured by the reduction of A Mcr . We note that 
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Fig. 6. (Color) Mean field ground-state phase diagrams for [Au(tmdt)2]. 
Parameters used are: fitted results for (a) room-temperature (T) structure, 
(b) low-T structure, and (c) room-T structure but with reduced Aj^ — > 
Aj, a — 0.1 eV. Dashed lines are guides for the eyes. The legend symbols are 
shown in (d); colors are appointed for distinctions of which orbitals showing 
spin moments as red (L), green (Ma), and blue (both L and Ma), while filled 
symbols represent insulating states. 



compared is about 0.05 eV while other parameters are only 
slightly changed, whereas the phase diagram bears a drastic 
change. 

The L-AFM phase stabilized in the large- Z7l, small- Um re- 
gion, is the remnant of the L-AFI phase in the room-T pa- 
rameters. Its ordering vector is the same and the amplitude 
of rrmcmetir moment is similar to that in the AFT nhase in 
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Fig. 7. (Color) Parameter dependences of expectation values of electron 
density, (rii), and spin density, |(n^) — for each orbital site in the 

lowest-energy MF solutions for Au(tmdt)2- (a)-(c) correspond to the cases in 
Fig. 6(a)-(c), respectively. The left (right) panel in (b) and (c) is for fixed ratio 
of £/m = Ul (Um = 2C/l). QWgpVIcr] denotes the magnetic ordering vector 
for the Mcr sites which changes at the parameters indicated by the broken 
lines. (g ma g[L] is always (tt : 0, 0), as seen in Fig. 6(d).) 



Fig. 6(a); For example, in the left panel of Figs. 7(b), it 
reaches 0.37/xb per site at Um = Ul = 0.8 eV. Neverthe- 
less, the system does not turn into an insulating phase, which 
is due to the orbital mixing: The occupation numbers are 
(nY a ) — 0.1 and (n\) ~ 1.45 for the low-T parameters, 
and (nf a ) ~ 0.2 and (n\) ~ 1.4 for the reduced A Ma case, 
which are noticeably shifted from the case of room-T param- 
eters. Then this L-AFM state in the large- Ul region can be 
considered as a 'doped dimer-Mott insulator' due to the mix- 
ing with Ma orbital. 

On the other hand, in the small- Ul, large- Um region, the 
Ma-AFM or Ma-FM state is stabilized; the L sites remain 
paramagnetic and then the system is naturally metallic. The 
magnitude of the moment appearing on the Ma sites is lim- 



ited by its occupation number, as seen in Figs. 7(b) and (c), 
about 0.1-0.2/iB- It is difficult to discuss the origin of each 
spin ordering vector for the Ma sites, g m ag[Ma], within our 
limited supercell size, since its small filling factor would typi- 
cally favor longer periodicity; In fact, it takes different values 
delicately depending on the parameters. However, we con- 
sider that the region where spin moments appear is reason- 
able. Comparing the phase diagrams in Figs. 6(b) and (c), the 
latter has larger region of phases with moments on the Ma or- 
bital, due to smaller A Mo . and therefore larger occupation in 
Ma. 

Even when we enter the region when both Um and Ul are 
large where spin ordering on the Ma and on the L sites coex- 
ist, the system still remains metallic. Here, the spin moments 
on the L sites can be large similarly to the case of the L-AFM 
phase. For example, in the left panel of Figs. 7(c), it reaches 
0.35-0.4/iB per site at Um = Ul = 0.8 eV. In this sense, this 
is also the remnant of the L-AFI phase in the room-T param- 
eters, stabilized in the similar interaction parameters. 

We note that in the case of the low-T parameters, there 
is a wide region where the ordering vector is g m a g [Ma] = 
(0,0,0), namely, ferromagnetic (Ma-FM) or ferrimagnetic 
state when the L sites also show spin moments (L&Ma-FM), 
is stabilized. This can be ascribed to the large density of state 
at ep, due to the lower edge of the ID band from the Ma or- 
bital. 16 In the room-T structure this edge situates just above 
ep, which becomes near it for the low-T parameters. The 
states with g m a g [Ma] = (0, 0, 0) also appears in the reduced 
A Ma case as well [see Figs. 6(c) and 7(c)]. 

3.2 [Cu(tmdt)2] 

The MF ground- state phase diagram for [Cu(tmdt)2] is 
shown in Fig. 8, while the parameter dependences of orbital 
occupations and magnetic moments are shown in Fig. 9. Al- 
though the Mtt orbital is off from the full occupation com- 
pared with the other cases (see Table I), we find no contribu- 
tion of it to the magnetic properties, therefore we can set them 
aside again. The occupations of the other two orbitals are now 
rather close to 1, i.e., 1/2-filling. Then, as discussed above, an 
AF insulating state with both Ma and L showing magnetic 
ordering (L&Ma-AFI) is seen to be stabilized in the region 
where both Ul and Um are large. Its spin pattern is the stag- 
gered one as was shown in Fig. 5(c). 

In the large- Ul, small- Um region, on the other hand, the 
L-AFM state is stabilized; vise versa, the Ma-AFM state is 
realized for small- Ul, large- Um', In these states the spin or- 
dering vector is the same as in the L&Ma-AFI phase [g mag = 
(7r, 7r, 0)]; these states are continuously connected. However, 
they cannot open an insulating gap since, although the AF or- 
dering provides a band splitting at ep in the magnetic orbital 
sector, the remaining orbital is paramagnetic with a Fermi sur- 
face; In other words, the spin ordering and band gap opening 
are orbital selective. The critical value for the magnetic states 
is smaller for the Ma-AFM, consistent with the fact that the 
Ma sites have smaller transfer integrals with ID structure. As 
shown in Fig. 9, the magnetic moment on each site can be 
large, due to the large filling factor near 1. In fact, the calcu- 
lations show a jump in the charge density on each site across 
the PM o AFM boundary (first order phase transition), which 
results in the occupation number for Ma and L to be closer to 
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Fig. 8. (Color) (a) Mean field ground-state phase diagram for [Cu(tmdt)2]; 
Parameters used are the fitted results in Table I. Dotted lines are guides for 
the eyes. The legend symbols are shown in (b); colors are appointed for dis- 
tinctions of which orbitals showing spin moments as red (L), green (Mcr), 
and blue (both L and Mcr), while filled symbols represent insulating states. 
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the PM phase. 

One point we should note is that, although the phase dia- 
gram is mostly dominated by the q mag = (7r, 7r, 0) state, we 
find many self-consistent solutions with very close MF ener- 
gies in the region where the Ma orbitals are magnetic. In these 
quasi-degenerate states, the spin ordering vectors are of the 
form q mag [Ma] = (7r, *, *) with * = or 7r, namely, only the 
2a periodicity is robust. This suggests that the spin exchange 
coupling between the moments on the Ma sites is essentially 
ID, which is consistent with the fitted results in Table I where 



the A bond along [100] has the largest Ma-Ma transfer inte- 
gral. This is also consistent with the first-principles calcula- 
tions comparing different magnetic orderings. 21 On the other 
hand, the L orbital always orders with g mag [L] = (7r, 7r, 0): the 
two-dimensional AF is stable in the MF level. 

4. Discussion 

In this section we compare our MF results with experiments 
on the magnetic orderings in this family, keeping in mind 
that in the calculations quantum fluctuations are neglected; 
additional effects of such fluctuations as well as the possi- 
ble strong correlation effect can be speculated on top of the 
MF results. Discussions on [Cu(tmdt)2] is followed by that 
on [Au(tmdt)2], since in the former we can have a consistent 
explanation of the experiments. As for the latter, we present 
several possible scenarios for the magnetic phase transition in 
this compound, in light of the considerations on its Cu analog. 

4.1 [Cu(tmdt) 2 ] 

Let us give an explanation for the experiments in 
[Cu(tmdt)2] based on our calculations, in association with the 
discussions in the literatures. 19-21 The experiments show an 
insulating behavior in the resistivity, and magnetic properties 
of a ID Heisenberg spin 5=1/2 system likely due to the 
pda -MO, namely, the Ma orbital. 

The insulating state in our MF calculations is realized for 
the L&Ma-AFI phase (Fig. 8) when both /7l and Um are 
large. We can consider this state as a multi-band Mott insu- 
lator where both L and Ma orbitals possess localized spins 
of effective 5 = 1/2 each. In the MF calculation this state is 
accompanied with a three-dimensional AF ordering. We can 
deduce the additional quantum effect as a spin singlet forma- 
tion in the L network: the rather strong dimerization on the 
B bonds can bring about the nonmagnetic ground state in the 
L subunit, as far as the other ti m are small enough. Then the 
active spin degree of freedom arises only in the Ma sector; 
They have a ID character as discussed, therefore consistent 
with the Mott insulating behavior with ID S = 1/2 chains. 
This situation is schematically shown in Fig. 10(a). 

Another scenario is that the system corresponds to the Ma- 
AFM phase in our MF phase diagram at the small- Ul and 
large- Um region (considering the relation Um > Ui), while 
other effect beyond our calculation brings about the insulat- 
ing behavior. In this MF solution, roughly speaking, the mag- 
netic ordering brings about a gap at ep for the Ma band, but 
the L bands remains metallic; This corresponds to the state 
found in the first-principles band calculation. 21 One possibil- 
ity is that, since, in this state, ep locates at the middle of the 
L bands where the band overlap is small, small perturbation 
might bring about a band gap: This is now a band insulator 
due to the dimerization in the L sector. Then the spin degree 
of freedom is only from the Ma orbital, again considered as a 
Mott insulator, showing the same magnetic behavior as above. 
The largest degree of dimerization in L sites among the cases 
listed in Table I is consistent with both pictures. 

4.2 [Au(tmdt) 2 ] 

There have been puzzling experimental data for 
[Au(tmdt)2], as mentioned in § 1. Our MF results for 
the room-T structure in § 3.1.1, which suggest only the L 
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Fig. 10. (Color online) Schematic density of states (DOS) for the possible 
situations for (a) [Cu(tmdt)2] and (b) [Au(tmdt)2]. The Mtt orbital is omitted. 
The band splittings compared to the situations in Fig. 4 are due to the forma- 
tion of Mott gaps, driven by Umo- for Mcr orbitals while for L orbitals by 
(a) ?7l and (b) the effective on-dimer Coulomb repulsion indicated by U^ f . 
Coherent peaks expected in such doped Mott insulating systems are omitted 
as well. 



sion as we discussed previously, 16-18 facing difficulties in 
explaining the experimental facts. Main problems were: (1) 
The AFM (spin-density-wave) state in the calculations shows 
small magnetic moment, whereas when is increased to 
achieve a large moment state the system enters the AFI phase, 
i.e., the dimer-Mott insulator; The experimentally observed 
metallic ground state with large magnetic moment could not 
be reconciled. (2) Discussions only involving the L sector are 
incompatible with the absence of sign of a phase transition 
in the resistivity at T A f- Spin-density- wave states due to the 
nesting of the Fermi surface should result in a change in 
transport properties. 

In clear contrast, on the other hand, the results of MF cal- 
culations in § 3.1.2, using the low-T parameters [Figs. 6(b) 
and 7(b)] as well as the reduced- values [Figs. 6(c) 
and 7(c)], show possible AFM states with larger magnetic 
moments. In these states, due to the orbital mixing between 
L and Ma, AF ordering cannot open a gap at e F , therefore 
the system remains metallic. Such an involvement of multi- 
orbital also brings about possible explanations for the absence 
of anomaly at Taf in the resistivity, suggesting that the mag- 
netic and transport properties are carried by different degrees 
of freedom. Below, let us propose two possibilities based on 
our results, considering the strong correlation effect in addi- 
tion. We ascribe the regions where magnetic moments arise in 
our calculations as doped Mott insulating states, as discussed 
in § 3.1.2. 

(i) {L: doped Mott insulator, Ma: PM}. This corresponds 
to the L-AFM state in our calculations. The L orbital forms a 
dimer-Mott insulating state but doped with holes, which are 
provided from the Ma orbital remaining in a PM state. The 
Mott gap is due to the effective on-dimer Coulomb interac- 
tion, as indicated by /7£ ff in Fig. 10(b) [case (i)]. 

(ii) {L: PM, Ma: doped Mott insulator}. In the MF calcu- 
lation, in the large-Z7M, small J7l region, Ma-AFM/FM states 
are stabilized. If similar picture for the Mott insulating nature 
in scenario (i) is applied, this results in a situation shown in 
Fig 10(b) [case (ii)]. The Mott gap formation in the Ma sector 
is adopted from the discussions on [Cu(tmdt)2]. 

In both cases, the doped Mott insulating character will gen- 
erate magnetism, while there exist carriers unchanged upon 
magnetic ordering. In the MF phases corresponding to case 



(i), large spin moments appear on the L sites, whereas its or- 
dering is expected to affect the transport properties since the L 
bands are mainly responsible for the conduction. On the other 
hand, case (ii) is favorable in the sense that the L bands are 
paramagnetic and carries the charge transport, and then the 
Ma sites are responsible for magnetism; such a picture has 
been proposed based on experimental considerations. 38 The 
weaker dimerization in the L sites than the Cu analog obtained 
in our estimated transfer integrals may be a factor to stabilize 
such metallic L system. However, in our calculation the value 
of the magnetic moment on the Ma is limited to be small 
by its occupation, i.e., the small electron doping. This does 
not agree with such an argument. Enhanced mixing between 
the two levels, approaching the situation in [Cu(tmdt)2] is ex- 
pected to bring about a possible reconciliation with the exper- 
imental facts, although within our estimation of parameters 
such a large mixing is not achieved. Future works including 
the strong correlation effect are needed for further investiga- 
tion of this possibility. 

Finally, the observed peculiarly high T A f is hard to discuss 
from our calculation only for the ground state, and evalua- 
tion of transition temperature at the MF level is usually not 
reliable to be compared with experiments, especially when 
strong correlation is involved. The low-dimensionality in both 
L and Ma orbitals when considering the intra-orbital transfer 
integrals is apparently incompatible with such high Taf, con- 
sidering the fact that the energy scale is similar with other 
molecular conductors showing lower transition temperatures 
for magnetic ordering in general. One possibility is that the 
doped carriers in mixed orbitals effectively make the spin- spin 
interaction large, especially even in the c direction where the 
original parameters ti m for among the L and Ma orbitals are 
small, making the system three-dimensionally coupled. Then 
the low-dimensionality embedded in each orbital can be re- 
leased to raise the critical temperature. 

5. Summary 

We have constructed effective models for single- 
component molecular conductors [M(tmdt)2] (M = Ni, 
Au, and Cu) showing a multi-orbital nature. The tight- 
binding parameters are obtained by fitting to first-principles 
band structures. The fragment molecular orbital picture 
leads us to a systematic view of this family: the interplay 
between a characteristic anisotropic electronic network and 
the orbital energy difference can tune the electronic states by 
different choice of M, especially that between the pda -type 
and p7r-type orbitals. By taking into account the Coulomb 
interaction, we discussed, for [Au(tmdt)2] and [Cu(tmdt)2], 
the mean-field phase diagrams and magnetic solutions of our 
effective model. In the former compound, we suggest that the 
mixing between the two orbitals can play a key role in resolv- 
ing their puzzling experimental results. On the other hand, 
in the latter, the existence of a multi-orbital Mott insulator 
is suggested which is consistent with the experimental facts. 
Both of these cases are distinctive examples in molecular 
systems. 
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